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The geometric analysis of the gyromotion for charged particles in a time-dependent magnetic field 
by J. Liu and H. Qin [Phys. Plasmas 18, 072505 (2011)] is reformulated in terms of the spatial angles 
that represent the instantaneous orientation of the magnetic field. This new formulation, which 
includes the equation of motion for the pitch angle, clarifies the decomposition of the gyroangle- 
averaged equation of motion for the gyrophase into its dynamic and geometric contributions. 



In a recent paper fl], Liu and Qin studied the dy- 
namics of the gyrophase of a charged particle moving 
in a time-dependent (but uniform) magnetic field B(t) = 
B x (t)x + B y (t) y + B z (t)z. Liu and Qin then obtained an 
expression for the equation of motion for the gyrophase 
in terms of the Cartesian components (B x , B y , B z ) and 
their time derivatives. By using this Cartesian formula- 
tion, Liu and Qin demonstrated the existence of an an- 
holonomic (i.e., path-dependent) geometric contribution 
to the gyrophase shift that depends on the evolution his- 
tory of the magnetic field during a gyro-period. 

In related previous work, Littlejohn 0, H| showed that 
the Hamiltonian theory associated with guiding-center 
motion in an inhomogeneous magnetic field [i[ displays 
an anholonomic geometric phase associated with the gy- 
rogauge vector R = Ve\ • e2 constructed from two basis 
unit vectors ei and % = b X % that span the plane that 
is locally perpendicular to the local magnetic unit vector 
b(t) = H(t)/B(t). Liu and Qin |l|, on the other hand, 
showed that anholonomic geometric phases also exist in 
the problem of the gyromotion of charged particles (i.e., 
not guiding-centers) in a time-dependent magnetic field. 

In the present Comment, we show how the analysis of 
Liu and Qin can be greatly simplified by using the polar 
angle 6 and the azimuthal angle p describing the orien- 
tation of the spatially-uniform magnetic field [3| , where 
cos 0(f) = B z (t)/B(t) and t&np(t) = B y (t)/B x (t). We 
also clarify the role played by the pitch angle in the evo- 
lution of the gyrophase. We now reformulate the analysis 
of Liu and Qin [1| concerning the case of time-dependent 
(but uniform) magnetic field B(f). 

First, we introduce the decomposition of the parti- 
cle velocity v in terms of its pitch angle A (i.e., v • b = 
|v| cos A) and gyroangle £ (i.e., dv/d( = v X b): 



cos A b + sin A c 



(1) 



where the unit vector c = a X b = da/d( depends ex- 
plicitly on £ and the speed v = |v| is a constant of the 
motion. 

Next, we obtain the equations of motion for A and £ 



by introducing the magnetic unit vectors [3] 

b = cos 9 z + sin 9 p 

6 = - sin z + cos p = db/89 

(p = — sin ip x + cos ip y = dp/ dtp 



(2) 



Here, the "radial" vector b points in the direction of 
the magnetic field at each point in space (for a time- 
dependent uniform field, this direction is only a function 
of time). Since the unit vectors satisfy b = 6 X <p, the 
unit vectors 9 = ei and (p = ?2 provide a possible choice 
of unit vectors (e^e-j) in the plane perpendicular to the 
unit magnetic vector b (i.e., tangent to the unit sphere). 



I. GYROGAUGE GEOMETRY 

With the choice ([2]) for the unit vectors (b,ei,e2), we 
construct the magnetic one-forms 

db-ei =60 ] 
db-e2 = sin# 6p (> (3) 
dei • % = cos 9 dp ) 

where d denotes an exterior derivative [H, [f| . Since the 
vectors (b,ei,e2) are orthogonal, we also have de; • b = 
— db-e; and dSj • = — dei-ej, for i, j = 1,2. These 
definitions are not unique, however, since a rotation of 
the perpendicular unit vectors (9,p) {9' ,(p') about 
the b-axis generated by the gyrogauge angle ip: 



cos ip sin ip 
(pi \ — sin-0 cos-0 
leads to the new one-forms 



(4) 



db-e^ = cos-0 d# + sin^ sini/; dp = u>2 dt 
db • eJ, = — sin-0 69 
cos 9 6<p -+ 



de'i •% 



it? sm'tp ap = uu'2 at 

sin6> cosi/j 6p = — w\6t 
6%p = lj 3 dt 



(5) 
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Equation §5§ shows the deep connection between the 
magnetic one-forms and the Eulerian angular frequencies 
(uji, u>2i W3), which are defined in terms of the Euler an- 
gles (tp, 9, tp— tt/2). We can thus imagine an infinitcsimaly 
thin symmetric top spinning about its axis of symmetry 
b = % (at constant 9 and tp with dtp/dt 7^ 0), which also 
undergoes precession (at constant 9 with dtp/dt 7^ 0) and 
nutation (at constant tp with d9/dt 7^ 0). 
If we define the gyrogauge one-form 

1Z = cos#dy, (6) 

and we denote its gyrogauge transformation as 

K' = K + dtp, (7) 

then the property of gyrogauge invariance 

dK' = dK (8) 

is guaranteed by the identity d 2 = (corresponding to 
the vector identity V X V = 0). The gyrogauge one-form 
© has a simple geometrical interpretation in terms of the 
solid-angle two-form [j| 



d9 A sin#d</? = — d^cos^dt^ 



dll. (9) 



Hence, according to Eq. (j8|), the solid-angle two- form (|9]) 
is a gyrogauge-invariant. In addition, using Stokes' The- 
orem [f|, the solid angle il = J v $ defined by the open 
surface V on the unit sphere is also expressed as 



n 



dTZ 



K. 



(10) 



I'D JdD 

where dV denotes the boundary of T>. 



where w c (t) = qB(t)/m denotes the time-dependent cy- 
clotron frequency and, using Eqs. pT|) - p^|) , we find 



db 
~di 



d9 . „ dtp ^ 

— 9 + sin 9 — tp 
dt dt Y 

/■ dd • /- • a dL P 
cosC — — sine smt> — 
s dt s dt 



* . dtp\ 

smc —r + cose smu — c, 
dt dt 



dc 
dt 



d( dtp 
— — cos 9 — 
dt dt 



sin C 



d9 
~dt 



cos ( sin 9 



dp 



(14) 



(15) 



We note that these expressions can also be expressed as 
db/dt — uji a — u>2 c and dc/dt = uj^a + u>2 b in terms 
of the Eulerian angular frequencies (wi, LJ2, W3) defined 
in Eq. with the substitution tp — > n/2 — £ (i.e., gy- 
romotion occurs as a counter-rotation about b). With 
these new definitions for the Eulerian frequencies, we 
easily recover the standard relations dui/dC = — 
dui2/d<^ = oji, and du^/dQ = 0. 

Equation (|13[) can be divided into two separate equa- 
tions of motion: 



dX 
17 



db 
~dl 

dc 

dt ' 



, db ^ 
cot A — • a. 

dt 



(16) 
(17) 



By using Eq. p^|) . Eq. ([TBI) yields the equation of motion 
for the pitch angle A: 



dX 
~dt 



sin£ 



dO 
~dt 



cosC sin 9 



dtp 
~dt 



(18) 



II. PITCH-ANGLE AND GYROANGLE 
DYNAMICS 



Next, we introduce the gyration unit vectors b = c X a 
in the plane perpendicular to b: 



a = cos £ 9 — sin ( tp, 
c = — sin ( 9 — cos ( tp. 



(11) 
(12) 



These definitions are gyrogauge-invariant (i.e., a' = a and 
c' = c) under the transformation £' = ^ + tp and Eq. (j4j . 
Using Eq. (UJ, the equation of motion dv/dt = w c vxb 
becomes 



• x~ x db dX^ 

= uj c sin A a + cos A — + — c 

\ dt dt , 



.' dc dX f- 
+ SmA l di ~ -dt h 



(13) 



While this equation is not considered by Liu and Qin 
[l|, it plays an important role in the evolution of the 
gyrophase [see Eq. (|2"0|) below]. We note that, since 
the angular velocities d9/dt and dtp/dt are gyroangle- 
independent (i.e., they represent the rate of change of 
the orientation of the magnetic field), the pitch-angle 
equation (fT5|) satisfies (dX/dt) = (which is valid for a 
uniform magnetic field), where (■ ■ • ) denotes a gyroangle- 
average. 

By using Eqs. (|H ]) -([I5 )) . on the other hand, Eq. (fT7|) 
yields the equation of motion for the gyroangle £: 



dt 



„ dtp 
+ cos 9 

dt 

J j-M . , . n dtp 
cot A cos C — — sin C sin 9 — 
V dt dt 



0J C 



cos 9 



dtp 
~dt 



uji cot A. 



(19) 



Equation (|19l) corresponds exactly to Eq. (7) of the paper 
1] by Liu and Qin, with the substitution £ — > — ( — tt/2 
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and 



B 



B X By - By B X 



B 2 



B 2 

v 



a d( ? 
— cos 9 ——. 

dt 



In Eq. (|19[) . uj c is described as the dynamical term by Liu 
and Qin, while cos 9 dtp/dt, which is clearly related to 
the gyrogauge one-form © , is described as the geometric 
term. The last term on the right side of Eq. (fl"9j) is de- 
scribed by Liu and Qin as an adiabatic term (because it is 
shown in the Appendix of Ref. [l[ to be one order higher 
than the geometric term when the magnetic field evolves 
slowly compared to the gyration period). Instead, we use 
the pitch-angle equation (fT8|) to write d(d\/ dt) / dC, = u>x 
(once again valid for a uniform magnetic field), so that 
Eq. p^|) is written as 



dX 
dt 



cos 9 



dp 
~dt 



d_ 



cot A 



dX 
~d~b 



(20) 



The third term in Eq. (1191) therefore appears as an exact 
gyroangle derivative, which disappears when Eq. (|19[) is 
gyroangle-averaged: 

' d(\ n dip 

dt) =UJc + C0S6 tt- 



(21) 



We note that Eq. (f2"Tj) is gyrogauge-invariant since, under 
a gyrogauge transformation generated by tp, we have 




d0' 
It 



d^ 
~dt'' 



dip 
~dt 



dip 
~dt 



which follows from the gyrogauge transformation ([7|). 
Upon gyroangle-averaging, the equation of motion (|21[) 
for the gyroangle is therefore decomposed in terms of dy- 
namical and geometric terms only. 



III. GEOMETRICAL CONTRIBUTIONS TO 
THE GYROPHASE 

We now follow Ref. [l| and use Eq. (|2~T|) to calculate the 
averaged gyrophase shift, denoted (AC), in one gyration 
period T = 2n/uj c (t): 

(AC) = dt' = (ACd) + <AC g )[C]. (22) 

The first term on the right side of Eq. (f2"2"j) denotes dy- 
namical gyrophase shift 

rt+T 

(ACd) = / LJ c (t')dt', (23) 



which equals 2n for a time-independent magnetic field. 
The second term on the right side of Eq. (l22j) . on the 
other hand, denotes the geometrical gyrophase shift 

<ACg)[C] = J f t+T cos 9(t')^dt' 



cos 9 dip = 71, (24) 
c Jc 



where the path C moves on the unit sphere from the 
initial point at 9(t) and <p(t), to the final point at 6(t + T) 
and ip(t + T). We note that the geometrical gyrophase 
shift (IM1) is path-dependent (i.e., it is anholonomic) since, 
by constructing the closed contour dV = C\ — C2 from two 
paths with identical end points, we find 



<AC g >[d] - <AC g )[C 2 ] 



TZ = / dTZ 

dv Jv 

• n + 0, 



(25) 



where ft denotes the solid angle enclosed by the open 
surface T> on the unit sphere. 

Lastly, we introduce the time-scale ordering on the evo- 
lution of the magnetic field 



d\nB 



dt 



T 



dp 
~dt 



T 



« 1, (26) 



where the magnitude and direction of the magnetic field 
are assumed to change on the same slow time scale com- 
pared to the gyro-period T. By inserting this ordering in 
Eqs. ([23"1) -([24" |) . we find the dynamical gyrophase shift 



(ACd) 



2tt 



T ■ 



dhxB 



dt 



(27) 



and the geometrical gyrophase shift 



(ACg)[C] 



dp 
~dt 



cos[9(t + r)]dr 



dt 



(28) 



which is just T times the instantaneous value of the ge- 
ometric term in Eq. (I2ip . and is of order e. Finally, Liu 
and Qin compute the average value of the adiabatic term 
in Eq. (|20p and show that it is of order e 2 . In our view, 
this contribution disappears upon gyroaveraging. 



IV. SUMMARY 

The formulation of the gyromotion in terms of the 
spatial angles relates the gyroangle dynamics with the 
motion of a spinning rigid body, through a natural 
appearance of the Euler frequencies. In addition, it 
emphasizes the role of the gyrogauge one-form TZ in the 
geometric interpretation of the gyromotion anholonomy. 
Lastly, the inclusion of the pitch-angle dynamics in our 
formulation shows that it is related to the adiabatic 
contribution to the gyroangle dynamics and it explains 
its adiabaticity. 

We conclude this Comment with a few remarks con- 
cerning a general magnetic field that is space-time- 
dependent (details will be presented elsewhere). In this 
case, Eq. (|2"Tj) is replaced with 



ds_ 1 dx 
T dt 2 dt 



(29) 
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where dei -?2 = cdt + R-dx Q. Here, two additional 
contributions appear. The first one involves the (Frenet- 
Serret) torsion of the magnetic-field line r = b • R = 
(dei/ds) -e 2 , with = ds/dt used in Eq. (J29j). The 
second one involves the twist of the magnetic-field lines 
T m = b'Vxb = d\/ds, defined as the rate of ro- 
tation (denoted by the angle x) °f a nearby field line 
about the magnetic field line represented by b [7, 8], with 



v\\ r m = dx/dt used in Eq. (|29|) . We note that the tor- 
sion and the magnetic twist can be comparable in some 
magnetic geometries Q. The torsion contributes to an 
anholonomic (path-dependent) geometric gyrophase shift 
(A£ T )[C] = j c rds while the magnetic twist contributes 
an holonomic (path-independent) geometric gyrophase 
shift (AC x )[C] = i/ c d X =\ Ax- 
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